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1 Introduction 


Two-dimensional quantum gravity has been an interesting laboratory for the study 
of fluctuating geometry. Many aspects have been understood by means of Liouville 
field theory, matrix models and the transfer matrix formulation of the theory. In 
particular, the transfer matrix formulation |l| has been useful for the analysis of what 
we will call “quantum geometry”, i.e. aspects of geometry which have no classical 
analogy. Surprisingly, such a situation appears already in pure two-dimensional 
quantum gravity. The partition function for pure two-dimensional quantum gravity 
where the volume of space-time is hxed to V is 

^V = /^gab]S(f^-V), (1) 

where [gab\ denote equivalence classes of metrics under reparametrization. With this 
dehnition the partition function for a hxed cosmological constant can be written as 
the Laplace transformation of Zy\ 

Za = JdV e-^^Zy. (2) 

From (|^ it follows that each geometry [gat] is assigned the same weight (one), i.e. 
there is no classical minimum around which it is natural to expand. This is why 
certain geometric aspects related to Zy (or Za) will be truly non-classical. A close 
analogy is found for the free relativistic particle. Let [P{x, y)] denote the equivalence 
class of paths from x G to y E R^, up to reparametrization invariance, and 
L{[P]) the length of the path in R^. The propagator of the free particle has a path 
integral representation closely analogous to (0)-(|l) 

GLix.y) = lv[P{x,y)] <5(L([P]) - L), GMix^y) = JdL e-^^^GL{x,y). (3) 

It is seen that each world-line contributes with weight one in the path-integral repre¬ 
sentation of Gl, precisely as each geometry did in the path-integral representation 
of Zy. It is well-known that a “typical” path [P{x,y)] has an (extrinsic) fractal 
dimension = 2. For instance, let us consider the ensemble of all equivalence 
classes of paths of length L. The corresponding partition function is 

Gl = Jdx Giix^y), 

and we can calculate 

{\x-y\)L = ^ JdxV[P{x,y)] 6{L{[P{x,y)]) - L) |a; - y] ~ (4) 

This is one of the simplest, but also most important quantum phenomena: as long 
as we address distances less than the inverse (renormalized) mass, it makes no sense 
to talk about any ordinary, one-dimensional path of the particle. Only for distances 
much larger than the renormalized mass, one can talk about an approximate classical 
path. 
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In the case of pure two-dimensional quantum gravity we have a somewhat sim¬ 
ilar situation, only will geometries and fractal dimensions refer entirely to intrinsic 
properties, with no reference to any embedding space. Let Sv{R) denote the average 
volume of a spherical shell of geodesic radius R in the ensemble of geometries with 
volume V. It can be shown that [@, §] 

+ (5) 


For any compact manifold of dimension d and a given, smooth geometry [gab] we 
have that 

SviR) ~ R'^-^ for R^O. (6) 

For an ensemble of geometries we call the power dh which appears instead of d for 
the intrinsic fractal dimension or the intrinsic Hausdorff dimension. From (j^) it 
follows that dh = ^ for pure gravity, rather than dh = 2 as one would naively have 
expected. Also, a calculation analogous to the one leading to (§) gives 

Wv ~ ( 7 ) 


which expresses that the average distance between two points in the ensemble of 
geometries only grows as and not as 

The precise definition of Sv{R) in pure quantum gravity is as follows: 


Sv{R) = /pm H jvs-v) I - R), (8) 

where ^ 2 ) denotes the geodesic distance between the points labelled by .^1 and 

,^ 2 - From the explicit calculation of Sv{R) in pure gravity we know that |^, ^ 

(9) 


where /(O) >0 and f{x) falls off like e for large x. Note that follows from 


From d^) it follows that Sv{R) can be viewed as a kind of reparametrization 
invariant two-point function between points separated by a geodesic distance R. 
The definition can be generalized to include matter fields. For a given metric Qab the 
reparametrization invariant partition function for matter will be denoted Z^[gab\, 
and it will appear as a weight in (|) and In this case it has not been possible 
to calculate dh by the same constructive arguments which led to = 4 for pure 
gravity. However, there exist arguments to be reviewed in the next section, 
based on the diffusion equation in quantum Liouville theory, which strongly suggest 
that dh{c) is a non-trivial function of c given by: 


4(c) 


^^/W^c+^/W^c 
v/25^+v/T^ ■ 


( 10 ) 


This formula agrees with the constructive approach for c = 0 and for c —> — cx) 
4(c) —>■ 2 as one would naively expect. As we increase c space-time becomes more 
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fractal until the analytic formula breaks down for c > 1. For c = —2 we have a very 
precise verihcation of (^) by numerical simulations However, for 0 < c < 1 the 
agreement with numerical simulations is not so good ilJ. We will return to this 
question below. 

While the dehnition of fractal dimension based on Sv{R) is in many ways natural, 
it is not the only one available. An alternative dehnition is based on diffusion and 
the dimension dehned in this way is called the spectral dimension. The dehnition 
has the advantage that it makes sense when dehned on “fractal structures” and we 
have just argued that a “generic” geometry in two-dimensional quantum gravity in 
a certain sense is fractal. For a hxed (smooth) metric Pab the dihusion equation has 
the form: 

( 11 ) 

where T is a hctitious dihusion time, is the Laplace operator corresponding to 
the metric gab and Kg{^, T) denotes the probability density of dihusion from ^ to 
,^0 in dihusion time T. If we consider dihusion with the initial condition 

eo; T = 0) = - eo) (12) 

V^(0 

it is well-known that Kg has the following asymptotic expansion for small T: 

e-n'=(€,&)/4r oo 

(13) 

r=0 

In particular the average return probability 

1 r I - 1 °° 

RP,(t) ^ ~ (14) 

where 

Ar = ^ /V^«r(e,0- 

The power rehects the dimension of the manifold, the heuristic explanation 
being that small T corresponds to small distances and for any given smooth metric 
short distances mean hat space-time. However, the dehnition is more general, and 
can be applied for dihusion in fractal structures, with the Laplacian appropriately 
dehned, as is well-known from the theory of percolation. From (|^ we have (for a 
smooth metric Qab) the classical result 

1 (D,(U<>)?K,(Uo-,T) ~T + 0 [r-), (15) 

irrespectively of d. 

Since the probability Kg is invariant under reparametrizations it makes sense to 
dehne the quantum average of Kg over all metrics: 

Kv{R]T) = Sv{R)V ^iJV9-V) Z^[9ab] 

X 6iDgi^,^o)-R) Kgi^,^o;T). (16) 
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By definition we have 


poo 

/ dRSv{R) Kv{R;T) = 1, 

Jo 

and fnrtherniore, the qnantnm gravity average of RPg{T) is 


RPv{T) = 




ab 


| 5 ( jvs - y) Zm W RP,{T) = Kr{0-, T). 


( 17 ) 


(18) 


It natnral to assnme that Kv{R]T) and RPy{T) have asymptotic expansions 
somewhat like (^) and (^). However, the powers of T which enter might be differ¬ 
ent from the canonical ones obtained for a hxed, smooth geometry. This sitnation 
is well-known from the stndy of diffnsion on hxed fractal strnctures (see [Q for a 
review). One operates with two different exponents. A dynamical exponent (or 
dimension) related to diffnsion (or random walk) on the fractal strnctnres and a 
strnctnral dimension, which we here identify with the intrinsic Hansdorff dimension]^ 
dh- The exponent 5^ is dehned by the mean-sqnare displacement after time T: 




'^‘2j6u 


(19) 


assnming that R{T) <C This means that the volnme covered by diffnsion 

after time T will be V{T) ~ {R{T))^ , and the probability that the random walk 
will retnrn to the origin shonld behave as: 


RPv{T) ~ 

Thns we have, by dehnition. 


+ o{T)) = ^^ 5 ^ 75(1 + o(T)). 


da 


‘J'dh 


( 20 ) 


( 21 ) 


Onr task in two-dimensional qnantnm gravity is to determine two of the three qnan- 
tities dh, dg and 5^. 

In the theory of diffnsion on fractal strnctnres it is nsnally assnmed (and well 
established nnmerically) that in the limit H—>^cx) K^{R-,T) has the following fnnc- 
tional form 

I ^ r R 


KooiR;T) = 






( 22 ) 


where Foo{x) falls off approximately as e“*“. Varions valnes of u has been considered, 
ranging from u = Ito u = 5^/ (5^ — 1). The fnnctional form (^) of conrse reprodnces 
(|I^) since 

(23) 


(fl”(r))„~ R" Kaa{R;T) ~ T”-"" 




®In the study of diffusion on fixed fractal structure one usually imagines the fractal structure 
embedded in . Thus one has an extrinsic fractal dimension Dh and intrinsic fractal dimension 
dh, the last one defined with respect to the “geodesic distance” of the fractal, which is defined from 
the shortest path between to points on the fractal. One usually has dh = i'Dh for some positive 
constant v. In the same way one has a relation similar to only with the distance Re(T) 

measured in R^, rather than intrinsically on the fractal: 

(i?|(T))^ 


The exponent Sw = f'Aw- 
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In the case of two-dimensional quantum gravity we want to consider a fixed volume 
V and average over all shapes. The original heat kernel expansion for a hxed geom¬ 
etry contains reference to powers of the curvature, but since we integrate over all 
geometries one expects that only reference to V will survive. We thus conjecture 


that 


VKv{R-,T) 


JL.] = JTf 

J'ds/2 Y2/dJ J'ds/2 Kylldh' Y2lds)' 


(24) 


where we have used (1^) to write 


F 


R 

T~ 1 



( ^ 

■ T ■ 

\VF 

-1 

_1 





(25) 


We expect the following boundary conditions on F and F: 


F{x,y) 


for ?/ —> cx). 

(26) 

F{x,y) 

0 

for y —>■ 00 , X > 0, 

(27) 

F{x,y) 

^ Fooix) 

for ?/ —> 0. 

(28) 


(^) results from the fact that Kv{R',T) 1/V for T ^ oo. We obtain (|^) from 
and (|?7|) because of (^). The above conditions are verihed in our numerical 
simulations. 

Finally the return probability for a hnite volume will be given by 


RPv{T) ^2/ds)‘ 

These are the scaling ansatze we will use in the following. Note that they imply the 
following: 


{R'^{T))y ~ for T^O, (30) 

(i?"(T))^ ~ for T^oo. (31) 


For any hxed, smooth geometry [gab] dh = dg = d, where d denotes the dimension 
of the underlying manifold (i.e. equal 2 in two-dimensional quantum gravity). After 
taking the functional average over geometries we know that dh changes, as already 
discussed. However, we will provide evidence that dg is unchanged and equal to two 
for all values c < 1 of the central charge c of the matter fields coupled to quantum 
gravity. In this context it is worth to recall that there exists a recent analytical 
argument in favour of this scenario 0: for Gaussian helds W^(,^i, (^ 2 ), /^ = 1, ■ ■ ■, T), 
coupled to two-dimensional gravity it is possible to derive the following relationship 
between the extrinsic Hausdorff dimension Dh of the surface X^((^i,^ 2 ) embedded 
in R^ and the spectral dimension dg, 


dg 


2Dh 
Dh + 2 


(32) 


One assumption going into this derivation is the scaling ansatz (p4D for i? = 0. 
Next, assuming that we can perform an analytic continuation of D from positive 
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integers to D G ] — cxo, 1[ one can appeal to Lionville theory and argne that Dh = oo 
for these valnes of D. Thns dg = 2 for this particnlar model. The nnmerical 
experiments reported in this article will provide evidence that the resnlt dg = 2 
is of larger generality than what was proven in |^, and also provide snpport for 
the scaling Note that for branched polymers (where Dh = 4) (|3^) resnlts in 

ds = 4/3 in agreement with . 

The rest of this article is organized as follows: In section we review shortly 
the diffnsion in quantnm Lionville theory and the derivation of (0), while section ^ 
presents the nnmerical evidence for scaling. Finally section ^ contains a discnssion 
of the resnlts obtained. 


2 Diffusion in Lionville theory 


Let *hn[fi'ab] be a general spin-less operator which depends only on the metric Qab^ is 
reparametrization invariant, and satisfies = A“”<hn[5'ab] at the classical level. 

The expectation valne of this operator in the context of two-dimensional qnantnm 
gravity conpled to a conformal field theory with central charge c is defined by 

{^n[gah])v = jT^bah] d( Jy/g-V) Z^lg^b] ^nbab]- (33) 

It follows from Lionville theory (see @] for details) that we have the following scal¬ 
ing: 


{■tnlSai.]),,. = {4„|S<.d>V . 


2ny/25 — c 

v/25^ + V25 - c - 24n■ 


(34) 


Consider now the diffnsion kernel Kg{^,^o;T) discnssed in the introdnction. The 
formal solntion is given as 




(35) 


We get the retnrn probability by setting ^ (after acting with and taking 

the average over all ^o- If we expand in T we obtain: 




(1 + TA, + --.) 


1 

7 ^ 


m - &) 


« 0 =? 


Let ns assnme the existence of a T' snch that 


(36) 


WRPxviT') = VRPviT). (37) 

From the assnmed scaling ansatz (0) it follows that 

rp! _ y2/dsrj^ _ yS-m/df^rj-^ (38) 

Since the scaling properties of the operator will change when dressed by two- 
dimensional quantum gravity, it is clear that one cannot maintain the combination 
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TAg in (P^D and (^) with T having it’s canonical dimension after averaging over all 
geometries. A better guess is obtained as follows: the average of the square of the 
geodesic distance travelled by diffusion at time T for a hxed geometry, as dehned in 
( 0 ). is again a reparametrization invariant object, and it makes sense to dehne the 
average in the ensemble of two-dimensional geometries weighted by Z^[gab]- Naively, 
one would expect 


/~T + o{r-), ( 39 ) 


the hrst term proportional to T coming from TAg if we use the expansion (PB|), and 


since we might expect a more general expression (|30D after averaging over geometries. 


it is natural to assume that one should consider the combination Then 


the hrst term in the expansion (|36[) would reproduce the behaviour (^OD, while (|36[) 
and (1^) and the scaling properties of Ag would allow us to conclude that 


dim 



dim 


y—a-ilai 


( 40 ) 


From (j^) and (|^ we arrive at the formula (|T^) : 


2ai ^a/25 — c -f i/49 — c 
0—1 — c 


( 41 ) 


However, note that this kind of argument does not allow us to determine the dimen¬ 
sion of T, i.e. ds or 5yj. This will be the purpose of the rest of the article. 


3 Numerical methods and results 



Figure 1: Finite size scaling of the return probability for the Ising model. 

We used dynamical triangulations in order to simulate conformal matter cou¬ 
pled to 2d quantum gravity on the lattice. In this approach, a triangulation T 
























Figure 2: Finite size scaling of the average distance travelled by the random walker 
for the Ising model. 


corresponds to an equivalence class of metrics [gab] in (0) and the volume V of 
spacetime is given by the number of triangles N in T. We used standard Monte 
Carlo techniques for unitary matter with c = 0 (pure gravity), 1/2 (Ising model) 
and 4/5 (3-states Potts model) and an effective recursive sampling technique for the 
(non-unitary) c = —2 model which constructs independent conhgurations. Details 
of the methods and models we used can be found in p, |^ . A certain number of con¬ 
hgurations was generated and the diffusion held Kjf{P, Pq;T), dehned on vertices, 
was evolved using the discretized version of (pd]) 


Kj,(P,P„-,T+l)^Yl 

j 


Kl{P,,P,-T), 


n{Pj 


iF/;(P,Po;0)=5p,Po, (42) 


where j runs over the neighbours of P and n{Pj) denotes the connectivity number 
of the vertex Pj. Then we obtain 

= . (43) 

dri^P, Pq) is the geodesic distance (shortest link path) between the points P and Pq 
and S')v(P) is the number of vertices at geodesic distance R from the point Pq. The 
return probability RPn{T) and the moments (P”(T))^ can easily be calculated 

RPn{T) = {kUPo,Po;T))^, (44) 

oo 

= Y.R''SN{R)K^iR;T) (45) 

R=0 

= (j2dr{P,PorKUP,Po;T)'j . (46) 
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c 

n = 0 

77, = 1 

n = 2 

dg 

b 

dh 

a 

dh 

a 

-2 

0 

1/2 

4/5 

2.00(3) 

1.991(6) 

1.989(5) 

1.991(5) 

2(5) 

4(5) 

4(4) 

5(5) 

3.58(13) 

4.09(23) 

4.08(32) 

3.99(24) 

0.6(3) 

1 .2(6) 

0.9(5) 

0.7(5) 

3.59(12) 

4.08(25) 

4.09(28) 

3.98(18) 

0.6(3) 

1.1(5) 

0.9(5) 

0.7(5) 



n = 

3 

n = 

4 

c 

dg 

a 

dh 

a 

-2 

3.55(9) 

0.5(3) 

3.53(8) 

0.4(2) 

0 

4.10(20) 

1.2(4) 

4.10(15) 

1.2(5) 

1/2 

4.10(25) 

1.0(4) 

4.11(23) 

1.0(4) 

4/5 

3.98(16) 

0.7(4) 

3.98(14) 

0.7(4) 


Table 1: The values for the spectral dimension dg and the fractal dimension dh 
obtained from calculating ^n{y) using hnite size scaling. The sizes of the lattices 
are = 2, 4, 8 and 16K triangles. 


N 

dg 

b 

dg{b = 0) 

128000 

1.980(14) 

1.0(7) 

1.9586(4) 

64000 

1.972(18) 

0.9(5) 

1.9400(5) 

32000 

1.958(32) 

0 .8(8) 

1.925(1) 

16000 

1.954(18) 

0.9(3) 

1.8949(7) 

8000 

1.938(34) 

0.8(4) 

1.865(2) 

4000 

1.934(58) 

0.9(4) 

1.826(3) 


Table 2: The values of the spectral dimension dg obtained from the small time 
scaling of <ho(l/) for the c = —2 model. 


It is more convenient to use (^Bf) for calculating (i?"(T))^. Evolving the held turns 
out to be an expensive procedure, so only one point Pq was chosen per conhguration. 
For this reason, the sampling of conhgurations in the case of unitary models was 
done sufficiently far apart so that they were essentially independent from each other. 

From the scaling relations one expects that RPn{T) and (P"'(T))^ 

will be functions of the scaling variables 


P T Cl 
Nl/dh ’ 

T + b 

N^/ds ’ 


(47) 

(48) 


where the “shifts” a and b are the lowest order hnite size corrections to scaling. 
The shift a has been used with great success in measuring correlation functions 
as functions of the geodesic distance B0.0 improving dramatically their scaling 
properties and making it possible to probe the fractal structure even for moderately 
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Figure 3: Small time behaviour of the return probability for the c = —2 model. 


small lattices. We expect such corrections to be necessary in our case as well. The 
scaling relations to be tested in the simulations are 


RPn{T) 

= ^^oiy), 

(49) 



(50) 


As we will see, from one obtains ~ 2 for all models. Using this value in 
( 0 ). one obtains values for dh consistent with the ones measured using different 
observables 10] verifying this way eq. (0. 

The functions ^n{y) for small y are expected to behave as: 

$o(l/) ~ (51) 

$n(l/) ~ n>0. (52) 

The validity of (^9D and (|50D was tested by collapsing the distributions for dif¬ 
ferent values of N ranging from 2-16K triangles. The method is described in full 
detail in §0- Measurement time grows as and this puts a severe limit on the 
maximum size of conhgurations possible to be studied. Measurements were per¬ 
formed on approximately 50000 configurations (10000-14000 for the 16K lattice). 
For the unitary models a configuration was obtained every 100 sweeps. One point 
Pq was randomly chosen on each conhguration. The best values for dg and d^ are 
recorded in Table [|. In the case of (i?"'(T))^, dg was hxed to be equal to 2, b was 
set to 0 and the fractal dimension dh as well the shift a were the free parameters to 
be tuned. The introduction of the shift a is crucial for these functions to collapse 
reasonably. In the case of RPhf{T) the only parameters involved are the T-shift b 
and the spectral dimension dg. It is worth mentioning that the collapse was done 
for a wide range of y (0.01-1) and that (x^/dof)min was considerably less than one 
in all cases (0.2-0.5 for approx. 45000 dof in each group). The errors quoted are for 
the range where y^/dof = 1. Figures ^ and show how well the scaling relations 
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Figure 4: Small time behaviour of the moments ^n{y) for the c = —2 model. N = 
4-128K from right to left for each distribution. The straight lines are oc 
ds = 2 and dh = 3.58. Slow convergence is observed as N ^ oo. 


hold in the case of the Ising model. Similar hgures can be obtained for the other 
models as well. 

Measuring dg and dh using the small time behaviour (^T]), (|5^ is more difficult. 
Larger lattices and more statistics are necessary for a sensible measurement to be 
made. For a detailed study, we conhned ourselves to the c = —2 model where it is 
easy to generate large conhgurations. Measurements were made on approximately 
80000 conhgurations (41000 for the 128K lattice). We evolved the diffusion held up 
to T = 1000 for lattices with N = 4iF-128iF. We checked that the results were 
consistent with the measurements we obtained from the unitary models, although 
with much less accuracy. 

In the case of $ 0 ( 2 /) the hts were performed by introducing the shift b and then 
by making a log-log plot for small y. The value of y^/dof was determined for a range 
of b from which we computed the best values of dg and b and their errors quoted in 
Table |l[ The small T cutoh Tmin = 7 was hxed for all volumes such that it would be 
the smallest Tmin giving x^/dof of order 1 for a reasonable range of T. The upper 
limit was hxed Uraax (be. T^a.x cc N). The hts were reasonably stable with diherent 
choices of T^in and Tmax- The values of dg for 6 = 0 for the same T-range are also 
shown in Table |l] for comparison. We see that b improves the value for dg quite a lot 
for the small lattices. In Figure we show graphically that (^) holds with dg = 2 
with very good accuracy. 

The above method is not so successful in the case of ^n{y) forn > 0. We observe 
large hnite size ehects entering in the calculation, which grow with n, as can be seen 
in Figure The straight lines correspond to the expected slopes and we see a very 
slow convergence as N ^ 00 . The hts, even for n = 0, do not yield stable values for 
dh and the results depend strongly on the range of T chosen. One has to throw away 
several small T points in order to obtain reasonable values of x^- Finite size ehects 
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enter in eq. (^B] ) throngh the assnmption that Sn{R) ~ (which we know that 

for the size of the snrfaces stndied is valid only for qnite small valnes of R) and from 
the assnmption that Kj\[{R;T) $ 0 ( 2 ;) where 2 ; = R/T^R^ which holds only for 

N —>■ 00 . 


4 Discussion 


The nnmerical resnlts reported above are two-fold: a corroboration of the conjectnre 
that dg = 2 for conformal matter coupled to two-dimensional quantum gravity, and 
a test of the scaling conjecture (pTf ) and The test of = 2 was two-fold. For 
c = —2, 0,1/2 and 4/5 a measurement of the return probability allowed a test of the 
functional form of VKv{0',T) in the form (p4D, and in this way a determination of 
dg- For a given central charge c it was done by “collapsing” the measurements for 
various V of the return probability as a function the single scaling variable 

y = ( 53 ) 

This is possible with impressive accuracy for a wide range of Id’s and T’s ii dg ^ 2 
as described above. The second, independent, test was only carried out for c = —2 
and concentrated on the small T dependence of Kv{Q]T). According to (P^) a fit 
to the power fall off should allow a determination of dg. This approach was used 
in the first systematic investigation of diffusion in the context of two-dimensional 
quantum gravity 0 . It does not allow a determination of dg with the same precision 
as the “collapse” method, but has the advantage, from the point of view of computer 
resources, that one only need to evolve the diffusion process for a small time interval. 
Again the result is ~ 2. 

The final test of the scaling form ( p^ is performed under the assumption that 
dg = 2. A measurement of the moments {R^{T))y allows a test of the i?-dependent 
part of the scaling hypothesis (P^). Again it is done by “collapse” of the measured 
distributions of {K^{T))y for various values of T and V and we find that their scaling 
is consistent with the existence of a scaling variable 


R 


X = 




(54) 


over all scales on the surface. This is in agreement with measurements on different 
correlation functions like the loop-length distribution function i 1,0 (from which 
one obtains e.g. Sv{R))- It is possible to perform such a collapse for a narrow range 
of dh- In this way one obtains an independent measurement of dh, compared to the 
one obtained in 0 . The agreement with the dh obtained by a direct measurement 
of the intrinsic Hausdorff dimension is perfect. Alternatively, the consistency of the 
results can be seen as a confirmation of (pT|). 

Summarizing, we have verified that with high accuracy dg = 2. Further, the 


scaling relation (p^ seems to be valid. Thus, we have a remarkable situation: a 
generic geometry which appears in the path integral in two-dimensional quantum 
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gravity, is fractal with an intrinsic Hansdorff dimension dh (which is a fnnction of the 
central charge c of the matter conpled to gravity). On snch a ensemble of geometries 
diffnsion is “anomalous”, i.e. 

+ for (55) 


rather that B?{T) ~ T as in ordinary diffusion on a fixed smooth geometry. However, 
this anomalous diffusion is counteracted by the fact that the geodesic distance itself 
has an anomalous dimension, and if the only measure of diffusion was the return 
probability, such a fractal space-time geometry would appear indistinguishable from 
an ordinary smooth two-dimensional space-time geometry. 

The values of dh measured by diffusion agrees with the values determined so 
far by direct geometric measurement i I, I, 0]- In particular one observes perfect 
agreement with (|^) for c = —2 and c = 0. However, for c = 1/2 and c = 4/5, 
i.e. in the case of unitary matter coupled to gravity, there is not a very impressive 
agreement. Thus we are still left with one of the few remaining puzzles in two- 
dimensional quantum gravity: is d/i = 4 for the central charge c G [0,1], or does it 
follow the prediction (|^) in this range of c (as seems to be the case for c < 0)? The 
fact that several independent ways of measuring dh agree, and fail to confirm (pTl), 
indicate that either (p!oD is not valid for c > 0, or there is a very general reason for 
the failure of the numerical simulations. One such reason could be that the volumes 
V considered so far are too small. Indeed, there have been arguments in favour of 
large finite size effects for c > 0 


Toll, but it is difficult for us to understand that one 


then should be able to measure critical exponents of, say, the Ising model coupled 
to quantum gravity with great accuracy, if this model suffers severe finite size effects 
for the same volumes when it comes to measurements of geometry. In particular, 
it is difficult to understand such a discrepancy between finite size effects on critical 
exponents and geometry when it is believed that it is the fractal geometry which is 
responsible for the change in critical exponents of the Ising model from the values 
in flat space to the KPZ values 
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